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Abstract 

We prove that a diagrammatic evaluation of the Kubo formula for the electronic transport conductivity due the exchange of 
bosonic excitations, in the usual conserving ladder approximation, yields a result consistent with the Boltzmann equation. 
In particular, we show that an uncontrolled approximation that has been used to solve the integral equation for the vertex 
function is unnecessary. An exact solution of the integral equation yields the same asymptotic low-temperature behavior 
as the approximate one, albeit with a different prefactor, and agrees with the temperature dependence of the Boltzmann 
solution. Examples considered are electron scattering from acoustic phonons, and from helimagnons in helimagnets. 
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1. Introduction 

Tlie earliest systematic treatment of electronic 
transport in solids was obtained by using the Boltz- 
mann equation. The latter is known to become exact 
in the limit of vanishing scatterer density. This implies 
that it becomes exact in the limit of low tempera- 
tures for scattering processes that involve excitations, 
such as phonons, that get frozen out in the limit of 
zero temperature. For most problems, the Boltzmann 
equation cannot be solved in closed form, but sys- 
tematic schemes of successive approximations exist; 
for instance, the method of variational solutions [1]. 
Usually, the lowest order variational solution yields 
the exact temperature dependence of the transport 
coefficients; one of the best known examples is the 
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Bloch-Griineisen law for the electrical resistivity p due 
to acoustic phonons, which states that the leading 
temperature dependence of p is given by p cx [1] . 

A different approach to transport phenomena uses 
the Kubo formula [2], which expresses the transport 
coefficients in terms of equilibrium time correlation 
functions that can be analyzed by means of many- 
body perturbation theory. For the electrical conduc- 
tivity, the relevant correlation function is the current- 
current correlation. It is known that an analysis that 
respects the conservation of charge requires a consis- 
tent treatment of self-energy and vertex-correction di- 
agrams. As a result, the transport relaxation rate is 
different from the single-particle relaxation rate, and 
the former is much harder to calculate than the lat- 
ter. The problem can be cast in the form of an inte- 
gral equation for the energy-dependent transport re- 
laxation rate l/rtr(e) that is not easy to solve [3]. An 
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approximation that ignores the energy dependence of 

l/rtr(e) [4] reduces the integral equation to an alge- 
braic equation and leads to a result that reproduces 
the Bloch-Griineisen result for the electrical transport 
rate due to electron-phonon scattering. This is in con- 
trast to the single-particle rate due to phonons, which 
is proportional to T^. However, because of the known 
energy dependence of the electronic self energy, which 
can be quite pronounced, this approximation has been 
questioned, and even the very validity of the law 
as derived from many-body perturbation theory has 
been called into question [4]. While such considerations 
are not consistent with what is known from the Boltz- 
mann equation, and thus should cast more doubt on 
the Kubo formula or its evaluation then on the Bloch- 
Griineisen law, they have led to a popular notion that 
the law is not firmly established. In any case, re- 
ducing the integral equation to an algebraic one is in- 
deed an uncontrolled approximation that can be justi- 
fied only by reference to the Boltzmann equation, not 
within the diagrammatic formalism itself. 

Scattering by excitations other than phonons leads 
to analogous questions. For instance, in the ordered 
phase of helical magnets such as MnSi the Goldstone 
modes associated with the helical order - helimagnons 
- lead to a single-particle relaxation rate proportional 
to T^''^, and the approximation mentioned above sug- 
gests a transport rate proportional to T''^^ [5]. While 
this result is extremely plausible by analogy to the 
electron-phonon problem, and can be checked by solv- 
ing the Boltzmann equation, it would nevertheless be 
very valuable to ascertain it directly within the many- 
body formalism. The purpose of this paper is to provide 
this missing link: Wc prove that the integral equation 
for the transport relaxation rate has a solution that is 
proportional to for the electron-phonon problem, 
and to T®''^ for the electron-helimagnon problem. 



2. Integral equation for the transport 
relcLxation rate 



ductivity a: 
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with rie the electron density and m-e the electron mass. 
If one defines 



(m = 0, 1), (2) 



then the energy-dependent single-particle relaxation 
rate To is given by 

ro(£) = -iVp £ du [n (^) + / (^)] V^'iu) , 

(3) 

with iVp the electronic density of states per spin at the 
Fermi surface, and n(x) and fix) the Bose and Fermi 
distribution functions, respectively. For a given Vq , 
Fo is an explicitly known function of T and e/T. The 
vertex part A is given as the solution of the integral 
equation 



A(e) = 1 - iVp 



du 



The transport relaxation rate rtr(£) = h.[e)/Vo[e) is 
the solution of 
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F 
where 

Ko,i{e,u) = -iV] 



du 



Ko{e,u) - Ki{e,u) 

Ftr(M) 



(5) 



4n^)+nr)J ro(e) 



(6) 

We will be interested in problems where the frequency 
in the spectrum ^"(p, u) scales as a power of the wave 
number, so that V{' (u) /Vo (u) — »• for m ^ 0. Since 
both the frequency u and the energy e scale as the 
temperature, this implies that, effectively, Ki/Kq — » 
for T — > 0. This observation will be important later. 



Consider electrons scattered by excitations that lead 
to an effective scattering potential with a spectrum 
V"{p, u), where p and u denote the wave vector and 
the frequency, respectively. Standard developments [4] 
lead to the following expression for the electrical con- 



3. Analysis of the integral equation 

From Sec. 2 we see that we need to analyze the in- 
tegral equation 
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(7a) 

where the first part of the kernel has the property 

/oo 
duKo{e,u) = l, (7b) 
-oo 

which follows from Eq. (3) . Equation (7a) is a Predholm 
equation of the second kind, and in what follows we 
will use the known properties of this class of integral 
equations [6] . We will refer to equations in Chapter 1 of 
Ref. [6] as (S xx), and to theorems by section number, 
(S Seen Theorem m). 

3.1. The homogeneous equation with kernel Kg 

We first consider the homogeneous integral equation 

/oo 
du Ko{e,u)ip{u). (8) 
-oo 

It follows from Eq. (7b) that A = Ao = 1 is an eigen- 
value of the kernel Kq with an eigenfunction ipo{e) = 1, 
(S 34,35). Note that it also follows from Eq. (7b) that 
f^^deduKo{e, u) = oo. This means that the sufficient 
condition (S 40ff) for an iterative solution of Eq. (7a) 
is badly violated, which suggests that iteration will not 
converge. This in turn suggests that the behavior of Ftr 
will be qualitatively different from that of Fq. 
It is also useful to consider the adjoint of Eq. (8), 



tp{e) = X / du Ko{u,e)ip{u). 



(9) 



From (S Sec.lO Theorem 9) it follows that for A = Ao = 
1 Eq. (9) has a nonvanishing solution; i.e., there exists 
a nonvanishing function ipo (e) that obeys 



•o{u) = r 

J — c 



deMe) Koie,u). (10) 



3.2. The inhomogeneous equation with kernel Kq 

We next consider the inhomogeneous equation 

ip{e) = r^+xf duKo{t,u)ip{u). (11) 

(S Socio Theorem 10) says that a necessary and suf- 
ficient condition for the existence of a solution of Eq. 
(11) is 

/.oo 1 

(12) 
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de —- ^{e) = 
lo(e) 



for every solution ■ip{e) of Eq. (9). An explicit check of 

this condition would require to explicitly find the solu- 
tions of Eq. (9), i.e., the left eigenfunctions of the ker- 
nel Kq, which is difficult. However, we will see in what 
follows that one can construct pcrturbative solutions 
of Eq. (7a) that do not exist as if i — > 0, which implies 
that Eq. (12) cannot hold for all solutions of Eq. (9). 

3.3. A perturbative eigenvalue of the full kernel 

Now we consider the following generalization of Eq. 

(7a), 

1 

ip{e) = z^rn +^ du[Ko{e,u) - aKi{e,u)] ip{u), 

J-0(£) 7-00 

(13) 

which reduces to Eq. (7a) for A = a = 1. Our goal 
is to use a as a small parameter. It follows from the 
observation at the end of Sec. 2 that working to leading 
order in a and then setting a — 1 will yield a result 
that is valid to leading order in T as T — > 0. 

We know from Sec. 3.1 that for a = 0, A = Ao = 1 is 
an eigenvalue of the kernel with eigenfunction ipo (e) = 
1. For the generalization of the eigenproblem of Sec. 
3.1 to the full kernel K = Kq — aKi, 



</?o(e) 



= Ao / 

J — c 



du K{e, u) ipo{u), 



this suggests a perturbative expansion 



Ao = A^°^+aA;^'+0(Q") 



(1) 



with Ao""* = 1, and 



(14) 



(15a) 



(15b) 



with 'fi'i^\e) = 1. Similarly, for the solution of the ad- 
joint eigenproblem, 

/oo 
de '4>o{e) K{e,u), (16a) 
-oo 

we write 

M^)=4'\^) + (^4'\^)+Oia^), (16b) 

where ipjy'\e) is the function that solves Eq. (10). Al- 
though we do not know this function explicitly, we have 
ascertained its existence in Sec. 3.1. 

Inserting these series into Eq. (13), multiplying with 
the function tpoie), and comparing coefficients yields 



Ao 



= l + a f 

J — c 

+0(a'). 



deVr(e)[Fi(e)/F, 



o(e)]/ r 

J — c 



de4°\e) 
(17a) 
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Here we have defined 

/oo 
duKi{e,u). (17b) 
-oo 

We see that the first order correction to the eigenvalue 

Ao is given by an average of the rate ratio Fi/ro with 
the zeroth order left cigenfunction Tp^'^ as the weight: 

Ao = 1 + Q (ri/ro)^(o) , (17c) 



A« = (ri/ro)^(o) , (i7d) 

where (/)„ = /^^ de /(e) V'(e)/ fZo ^'(e)- 

3.4. A perturbative solution of the full integral equation 

The desired solution of Eq. (13) can now be written 
in terms of a resolvent R{e, it), see (S 46), 

1 f°° 1 

= r77^ + / '^^ r-HA ' ^^^^ 

io(e) J-oo '■oW 

and the resolvent can be expanded in a Laurent series 

(S Sec.8), 

R{e,u) = ^^^^+O{{l-Xof). (19) 

Since 1 — Ao = 0{a), keeping only the leading term in 
the Laurent series is consistent with working to lead- 
ing order in a. The residue a_i satisfies the integral 
equation (S Sec.8) 



o-i(£,m) = Ao / dt K{e,t) a-i{t,u) 



-r 

J —C 



dt Ko{e,t) a-i{t,u) + 0{a). (20) 



If we define 

/oo 
du a_i(€, u)/ro{u) (21a) 
-oo 

we have 

dt Koie,t) A-i{t) + 0{a). (21b) 

-oo 

To lowest order in a, this is just Eq. (8) with A = 1, 
which has a solution 

A-i{e) = const. + 0(a) = A-i + 0{a). (21c) 

The solution of the integral equation can now be writ- 



ten 



^(e) = — 4t + 0("")- 



(22) 



We see that to leading order in the Laurent expansion 
in Q, (fi does not depend on e. We also see that the 



solution does not exist for a = 0, as we had anticipated 
in Sec. 3.2 above. 

The remaining task is to determine the constant 
A-i, which is not fixed by the integral equation (21b). 
Let us consider the functions <^(e) and l/ro(e) as el- 
ements \ip) and |l/ro), respectively, of a linear space, 
and the kernel K a linear operator in that space. The 
integral equation (13) then reads 



(23) 



Multiplying from the left with the left eigenfunction 
{ipo\ of K, we have 



{Vo|¥') = (^/'o|l/ro> + (Vo|^|¥'> 
= {i/>o|l/ro> + Ao(Vo|<^>. 



(24) 



From Eq. (22) we know that, to lowest order in a, \ip) 
is a multiple of Ii^q"^)' 

l'^> = ^ 4f l^r> + 0(a°). (25) 

To leading order in a we thus have 

(1 - Ao) (Vri<^r> [1 + 0(«)] = (4°'|l/ro>, (26) 
which yields 



A-i = 



(27) 



Using the notation defined in Eqs. (17) we have 

1 (1/ro) ,(0) 

^ie) = ^——^+0{a"). (28) 

Finally, using the observation made after Eq. (13) we 
can put a = 1 and write 



vie) 



<i/ro>,(o) 

^+0(r;). (29) 



<ri/ro)^(o, 



This is the desired solution of the integral equation (7a) 
to leading order in the function ri(£) defined in Eq. 
(17b). The only assumption we have made is that the 
left eigenfunction tp^'^ (e), which we have not explicitly 
determined, falls off sufficiently fast for large e to ensure 
the existence of the averages in Eq. (27). 

The above result is exact. Let us compare it with 
the popular approximation [4] that replaces ip{u) under 
the integral in Eq. (7a) by <^(e). This turns the integral 
equation into an algebraic one, and the solution is 



^(e)«l/ri(e). 



(30) 
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This has the same structure in terms of Fi as the exact 

solution (28), and leads to the same temperature de- 
pendence of the conductivity, as we will see. However, 
it is uncontrolled and cannot be justified other than by 
a comparison with the exact solution provided above, 
or with a solution of the Boltzmann equation. 



4. The electrical conductivity 

4.1. The Drude formula 

We now can write the conductivity as follows. From 
Eq. (1) we have 

-2 de 1 



1 



(31) 



2m, J _^ 4T cosh2(e/2T) Ftr(e) 

Equation (29) shows that, to lowest order in Fi, Ftr is 
independent of e and given by 

rtr= +0{V\). (32) 

We thus have as our final result the Drude formula 



rieC Ttr 



with a transport relaxation rate 
1 



— = 2Ftr . 

Ttr 



(33a) 



(33b) 



Here Ftr is given by Eq. (32). It depends on Fo from Eq. 
(3), Fi from Eq. (17b), the kernels Kq and K\ defined 
by Eq. (6), and the left eigenfunction tp'^^^ of Kq that 
is defined as the solution of the equation 

/oo 
deVr(6)i^o(e,M). (33c) 
-oo 



where Co and Ci are positive constants. For acoustic 
phonons, the exponents rim are [7] 

no = 3, ni = 5 (phonons). (35a) 

For helimagnons, the corresponding values axe [5] 

710=3/2, m = 5/2 (helimagnons). (35b) 
For the relaxation rates Fo,i this yields 

F^ (e) = T"- 7^ (e/T) , (m = 0, 1) , (36a) 

where 



7r 



/•OO 

,{x) = Cm du w"-"-' [2n(w) + f{u + x) 
Jo 

+f{u-x)]. (36b) 



Finally, it follows from Eqs. (6) and (33c) that V'o"^ (c) 
is a function of e/T only. For the transport relaxation 
time, and hence for the conductivity, we thus have the 
a temperature dependence from electron-phonon scat- 
tering: 

IT (X Ttr oc 1/T^ (phonons), (37a) 

This is the well-known Bloch-Griineisen law [1]. The 
same power law is obtained by means of the uncon- 
trolled approximation (30) [4] , but the prcfactor is dif- 
ferent. For scattering from licliniagnoiis wc find 



(T oc Ttr tx l/y'/"^ (helimagnons). 



(37b) 



Again, the uncontrolled approximation yields the same 
l/ys/a temperature behavior [5]. The above consid- 
erations ascertain that this is indeed the exact low- 
temperature contribution to the electrical conductivity 
due to electron-helimagnon scattering. 



5. Conclusion 



4.2. Two physical examples 

We finally illustrate our result with two examples. 
One is the well-known result for the conductivity due to 
electron-phonon scattering, the other, the conductivity 
due to the scattering of helical magnon excitations in 
helimagnets. In both cases, the effective potentials V'o"! 
that determine the kernels have the following power- 
law form in the limit of low frequencies: 

VZ{u) = ^ luf"^-' sgnu, (m = 0,l), (34) 



In conclusion, we have considered the Kubo formula 
for the electrical conductivity due to the scattering 
of electrons by bosonic excitations. We have given an 
exact solution of the integral equation for the vertex 
function that determines the transport relaxation rate. 
This exact solution has the qualitatively same behav- 
ior as the result of a popular, but uncontrolled, ap- 
proximation. For the scattering of electrons by acous- 
tic phonons this yields the usual Bloch-Griineisen law 
for the electrical resistivity, p oc T®, and for the he- 
limagnon contribution to the conductivity in helical 
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magnets it confirms the result p oc T^/^ that had been 

previously obtained by means of solving the integral 
equation approximately. All of these results are consis- 
tent with what is known from solutions of the Boltz- 
mann equation. 
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